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We propose a new design for a quantum bit with four superconducting islands in the topology 
of a symmetric tetrahedron, uniformly frustrated with one-half flux-quantum per loop and one-half 
Cooper-pair per island. This structure emulates a noise-resistant spin-1/2 system in a vanishing 
magnetic field. The tetrahedral quantum bit combines a number of advances such as a doubly 
degeneracy ground state minimizing decoherence via phonon radiation, a weak quadratic sensitivity 
to electric and magnetic noise, relieved constraints on the junction fabrication, a large freedom in 
manipulation, and attractive measurement schemes. The simultaneous appearance of a degenerate 
ground state and a weak noise sensitivity are consequences of the tetrahedral symmetry, while 
enhanced quantum fluctuations derive from the special magnetic frustration. We determine the 
spectral properties of the tetrahedral structure within a semiclassical analysis and confirm the results 
numerically. We show how proper tuning of the charge-frustration selects a doubly degenerate 
ground state and discuss the qubit's manipulation through capacitive and inductive coupling to 
external bias sources. The complete readout of the spin-components o~i, i = x,y,z, is achieved 
through coupling of the internal qubit currents to external junctions driven close to criticality during 
the measurement. 



I. INTRODUCTION 

Superconducting solid-state qubits (short for quantum 
bits) are promising candidates for the future construc- 
tion of quantum information processors. They appear in 
a variety of designs: in the charge-versioni2^ the quan- 
tum information is stored in the number of excess Cooper 
pairs residing on a small superconducting island — this 
design requires fabrication of ultra-small structures and 
is susceptible to charge noise. In the flux- 4 and phase-^i 
versions the information is encoded in the current state of 
the device — this is a macroscopic variable susceptible to 
flux noise. The new design by Vion et al£ is 'in between', 
with the energy scales for the charge- (Ec = e 2 /2C) and 
the phase- (Ej = & I c /2ttc) degrees of freedom roughly 
balancing one another (here, C and I c are the capacitance 
and the maximal current of the device, $o = hc/2e is the 
flux quantum) ; its ground state is non-degenerate and the 
limitations are close to those of the charge device. 

The novel tetrahedral qubit design we propose be- 
low operates in the phase-dominated regime and exhibits 
two remarkable physical properties: first, its non-Abelian 
symmetry group (the tetrahedral group Td) leads to the 
natural appearance of degenerate states and appropriate 
tuning of parameters provides us with a doubly degener- 
ate groundstate. Our tetrahedral qubit then emulates a 
spin-1/2 system in a vanishing magnetic field, the ideal 
starting point for the construction of a qubit. Manipula- 
tion of the tetrahedral qubit through external bias signals 
translates into application of magnetic fields on the spin; 
the application of the bias to different elements of the 
tetrahedral qubit corresponds to rotated operations in 
spin space. Furthermore, geometric quantum computa- 
tion via Berry phases^ might be implemented through 
adiabatic change of external variables. Going one step 



further, one may hope to make use of this type of sys- 
tems in the future physical realization of non-Abelian 
anyons, thereby aiming at a new generation of topological 
devicesifliiiiiSii 3 . which keep their protection even during 
operation^. 




FIG. 1: (a) Tetrahedral superconducting qubit involving four 
islands and six junctions (with Josephson coupling Ej and 
charging energy Ec); all islands and junctions are assumed 
to be equal and arranged in a symmetric way. The islands are 
attributed phases (pi, i — 0, ... ,3. The qubit is manipulated 
via bias voltages Vi and bias currents ij. In order to measure 
the qubit's state it is convenient to invert the tetrahedron 
as shown in (b) — we refer to this version as the 'connected' 
tetrahedron with the inner dark-grey island in (a) transformed 
into the outer ring in (b). The measurement involves ad- 
ditional measurement junctions with couplings E m 3> Ej 
on the outer ring which are driven by external currents 7 m 
(schematic, see Fig. 6 for details); the large coupling _B m ef- 
fectively binds the ring segments into one island. 

The second property we wish to exploit is geomet- 
ric frustration: In our tetrahedral qubit discussed below 
it appears in an extreme way by rendering the classi- 
cal minimal states continuously degenerate along a line 
in parameter space. Semi-classical states then appear 
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only through a fluctuation-induced potential, reminis- 
cent of the Casimir effect and the concept of inducing 
'order from disorder ilfitli. The quantum-tunneling be- 
tween these semi-classical states defines the operational 
energy scale of the qubit, which turns out to be unusu- 
ally large due to the weakness of the fluctuation-induced 
potential. Hence the geometric frustration present in our 
tetrahedral qubit provides a natural boost for the quan- 
tum fluctuations without the stringent requirements on 
the smallness of the junction capacitances, thus avoid- 
ing the disadvantages of both the charge- and the phase- 
device: The larger junctions reduce the demands on the 
fabrication process and the susceptibility to charge noise 
and mesoscopic effects, while the large operational en- 
ergy scale due to the soft fluctuation-induced potential 
reduces the effects of flux noise. Both types of electro- 
magnetic noise, charge- and flux noise, appear only in 
second order (cf. also Ref. 0). 

In the following (Sec. 2), we first introduce the struc- 
ture of our tetrahedral qubit and then find the low-energy 
part of its spectrum. We proceed in three steps, begin- 
ning with the (highly degenerate) classical solution; sub- 
sequently, we demonstrate how the fluctuation induced 
potential reduces this degeneracy to three semi-classical 
states and finally, we analyze the tunneling between them 
in order to arrive at the final answer for the phase dom- 
inated regime. We confirm this solution with the help of 
numerical calculations and extend it to the charge dom- 
inated regime. The inclusion of external fields breaking 
the (tetrahedral) symmetry of the device prepares the 
discussion of the qubit 's manipulation schemes (section 
3). In section 4, we discuss various measurement schemes 
and end with the conclusions in Sec. 5. A first account 
on part of this work has been given in Ref. Il8l 



II. TETRAHEDRON 

A. Device structure 

Consider the planar structure made from four super- 
conducting islands interconnected via six (conventional) 
Josephson junctions with equal couplings Ej\ connecting 
each island with all the others produces the topology of 
a tetrahedron, see Fig. 1(a). The islands are numbered 
through '0' to '3', with the island '0' residing in the cen- 
ter, and are assigned phases 4n, i — 0, . . .3. The three 
triangular loops are small (i.e., Ej -C ($o/27t) 2 /£a with 
La the loop inductance), allowing us to neglect fluxes 
induced by currents flowing in the structure. In the ab- 
sence of external magnetic fields, the classical energy of 
this arrangement is given by the sum Vq = J2i<j Ej [1 — 
cos(0jj)] with the difference variables 0,j — 4>j ~ A 
slightly modified version of this device with the inner 
island converted into an enclosing ring is shown in Fig. 
1(b) — choosing appropriate parameters, this variant ex- 
hibits the same physical properties as the original isolated 
tetrahedral qubit. In addition, this second design ideally 



lends itself for measurement of the qubit state. Below, we 
treat the ring as one connected island; the strong junc- 
tions with E m S> Ej used in the measurement process 
will be discussed later. 

Next, we bias the structure through an external mag- 
netic field, frustrating each sub-loop with a flux $o/2 
(these are three triangular sub-loops in the isolated ver- 
sion of Fig. 1(a) and 4 sub- loops in the connected ver- 
sion of Fig. 1(b)). We include the effect of this flux in 
a symmetric gauge by adding the phase n to each of the 
difference variables 0y-; the energy (up to a trivial con- 
stant; we measure phases with respect to the phase 0o of 
the central/ring island) 

Vx — Ej [COS 4>\ + COS 4>2 + COS 03 (1) 
+ COs(0 3 - 2 ) + COs(0i - 03 ) + COs(0 2 - 01 )] 

then is minimized (with V„ = ~2Ej) along the lines 

03 = ±7T; 02 - 01 = ±7T; 1p3 = 01 + 02 € [-7T, 7r] (2) 

and their analogs obtained by cyclic permutation 3 — > 
1 — > 2 — > 3. These minimal-energy lines run in the planes 
of the cube [— 7r,7r] 3 defined in phase space {01,02,0a}, 
see Fig. 2(a). The huge (linear) classical degeneracy can 
be easily understood via a reformulation of the poten- 
tial m in terms of the complex variables Zk = exp(i0fc), 
2V V = Ej[\ X)fc=o z k\ 2 ~ 4]; this expression is minimal for 
^ k Zk = 0. The two conditions defined by this (com- 
plex) equation imply that one out of the three variables 
0fc can be freely chosen, thus defining lines of minimal 
potential energy. Of particular relevance are the min- 
imal energy configurations O;, i — 1,2,3 on the cube 
edges where two minimal-energy lines join. These con- 
figurations involve two opposite junctions with a phase 
difference 4>ij — 0, while the remaining junctions involve 
maximal phase differences 7r; following three consecu- 
tive segments on the cube, these minimal states rotate 
through 27r, see Fig. 2(a) (here, we have included the 
bias phases 7r on each link; in the original variables these 
states involve 2 strained junctions with phase difference 
7r and 4 unstrained ones with phase difference 0; note the 
absence of currents in these configurations). 

We account for the quantum dynamics in the array via 
the capacitive term 



T = 



4e 2 



E 



o 2 + ^ 



Co, 
2 



in the Lagrangian C = T — V^; here, 



(3) 

Cj denotes the 
and Co are the 



capacitance of the junctions, while C g 
capacitances to the ground of the islands i = 1,2,3 and 
of the center/ring island. An additional term ft^ift^i 
appears when charges 2eqi are induced on the islands; we 
will discuss the effect of this topological term later. For 
the isolated tetrahedron of Fig. 1(a) we have Co = C g ; 
going over to center of mass (4$ = 0j) and relative co- 
ordinates (0j = 4>i — <!>) the difference variables 0^ = 0y 
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FIG. 2: (a) Continuous degeneracy of the classical minimal 
energy states. We identify 3 highly symmetric minimal states 
Oi, i = 1,2,3 (white, black, and grey dots) involving two 
opposite junctions with a phase difference <j>ij — and four 
remaining junctions with (/>y = tt. Three families of lines 
(solid, dashed, dotted) connect these states around the cube 
[— 7r, 7r] 3 where the energy is continuously degenerate with a 
value V-n — —2E.j. The (non-orthogonal) coordinates \i are 
directed along the minima, e.g., xi = 0, X2 = — tt, X3 £ [0, 7r] 
along the reference line 7 along which the change in the po- 
tential Vn is shown for the values X3 = 0,7r/2,7r. Quantum 
fluctuations induce potential barriers along the minimal lines, 
separating the classical states Oi and transforming them into 
semi-classical states \Oi). (b) Tunneling between the states 
\Oi) establishes the low-energy spectrum of the tetrahedron. 
The tunneling process connecting the states \Oi) and \0%) 
involves two trajectories 7 and 7' where a fluxon $o/2 cuts 
through 4 junctions, flipping all of them by -k. The phase 
difference between the two trajectories (the solid line corre- 
sponds to 7) produces the Aharonov-Bohm-Casher phase ob- 
tained when taking the fluxon $0 around the islands '1' and 
'2'. 



B. Semi-classical Analysis 

Next, we account for quantum fluctuations associated 
with the degenerate line-minima. It is convenient to 
introduce the (non-orthogonal) variables Xk — (4>i + 
4>j ~ 0fe)/2 with i,j,k £ {1,2,3} and mutually differ- 
ent, where both the kinetic and potential energy terms 
aquire a simpler form, T re i = (h 2 /4E c )[xl + X2 + X3] an d 

K = 2Ej [COS Xl COS X2 + COSX1COSX3 + COSX2COSX3]- 

The new coordinates are directed along the potential 
minima which are parametrized by fixing two coordinates 
to and ±7r and have the third run through the interval 
[0, =ftt]; e.g., the reference line (denoted by 7) connect- 
ing O2 with 0\ in Fig. 2(a) is parametrized by xi = 0, 
X2 = — tt, X3 € [0j7r]. I n the vicinity of 7 the potential 
energy takes the form V„ ~ Ej[— 2 + Sx 2 (l — cos X3) + 
5x1(1 + COSX3)]; the two fast oscillatory modes xi an d 
X2 appear with a curvature which depends on the adia- 
batic coordinate x = X3> see Fig- 2(a). Their zero-point 
fluctuations produce an induced potential 



Vf(x) = ^NW + ^W], 

with the frequencies of the fast modes 



wi :2 (x) = w fv /[l±cos(x)]/2 



(4) 



(5) 



and the frequency scale ujf = \[EEjEc / h. Near the edges 
X = 0, 7r one of the modes goes to zero, as is clear from the 
potential shape shown in Fig. 2(a), and we have to refine 
our analysis. We then remain with only one fast and 
two slow modes. Expanding V„ around the point O2 = 
(0, — 7r, 0) (in Xi coordinates) we arrive at the potential 
V* = (Ej/2)[-4 + ASxl + SxiSxl ~ 6x1(6x1 + Sxl)}; 
integration over the fast mode X2 and transformation to 
momenta provides us with the Hamiltonian 



Hi 



-E c 



dxl 



dxl 



E 



6x\5xl ~ k(6xI + 6x1) 



(6) 

where k = (6xi) = (Ec /^Ej) 1 ' 2 . Dimensional analysis 
tells that the low lying levels of this quartic anisotropic 
oscillator are of the order of 



2 , with the capacitive en- 
The identical 



pick up an additional capacitance Co/4 and the corre- 
sponding part of the kinetic energy can be written in the 

fonnT IBl = (fi 2 /16£o)Ei<y'M 
ergy E c = e 2 /2C and C = Cj + C /4 
expression for the kinetic energy is obtained for the in- 
verted tetrahedron in Fig. 1(b) if we choose a large self- 
capacitance Cq/Cj — » 00 for the ring and a small one 
Cg/Cj — > for the other islands, then E c = e 2 /2Cj; 
this limit describes the inverted tetrahedron connected 
to the outside world via (large) superconducting wires 
fixing the phase (f>o — 0. The following discussion applies 
to both designs of Fig. 1; we first assume that Ec <C Ej, 
placing the array into the phase-dominated regime. 



n = u; i [ — 



1/6 



< UJf, 



(7) 



the numerical factors determining the exact positions 
of the groundstate and of the non-equidistant higher 
levels have to be determined numerically and the re- 
sults are summarized in Table [I] the groundstate en- 
ergy is accurately described by the expression Oo/O « 
0.311 - 0.129(E C /Ej) 1 / 6 , where the scale dependence of 
the correction easily follows from first order perturbation 
theory in the term —Ejk(6x\ + ^xi)/2- 

Repeating this analysis for the other classical line- 
minima, we arrive at three distinct quantum states \Oi) 
(at energies —2Ej + huii/2 + MIq) associated with the 
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TABLE I: Ground- and excited states energies near the min- 
ima Oi. The first excited state at Q 1 is doubly degenerate. 
The parameter v = u^(Ej/Ec) 1 ^ 3 quantifies the susceptibil- 
ity of the semi-classical ground states \Oi) to applied fluxes, 
cf. l|23|l . 



energy of the two fast oscillating modes, 



Ej/Ec 


Q.o/0. 


Qi/O 


n 2 /Q. 


v * 


oo 


0.311 


0.70 


0.99 


1.00 


1000 


0.271 


0.57 


0.69 


1.23 


100 


0.250 


0.51 


0.58 


1.53 


10 


0.225 


0.40 


0.43 


2.50 



three classical zero-current states Oi described above. 
These isolated quantum states are generated through 
a fluctuation-induced potential reminding about the 
Casimir force between metallic plates^ or the van der 
Waals interaction between neutral atomsii This is just 
the mechanism producing order from disorder originally 
proposed by Villaini£*i£ where the huge classical ground 
state degeneracy (which does not follow from the sym- 
metry properties of the system) is removed by quantum 
fluctuations, here, selecting the three points \Oi) as new 
ground states. 

The low-energy spectrum near the points Oi exhibits 
non-equidistant levels f2i even deep in the semiclassical 
regime, allowing for the use of the tetrahedral struc- 
ture as a simple J osephs on junction qubit of the type 
introduced in Refs. 20 21. Moreover, both symmetry ar- 
guments and the numerical data tell that the first ex- 
cited level is twofold degenerate, such that we effectively 
deal with a spin 1 system with an easy-plane anisotropy 
H s =i = (fix-noJSj. 

Before going to the full quantum description, let us dis- 
cuss the above semiclassical version of the device, as it 
exhibits a number of interesting features by itself. First, 
the potential J3J defines a doubly-periodic junction^ 
with two distinct minima. The potential Vf(x) can be 
mapped out experimentally through the measurement of 
the Josephson current Ij (x) that can be pushed through 
the structure. E.g., fixing the phase fa between the 
central island and the island '2' via a flux-biased ex- 
ternal loop defines the two classical minimal solutions 
(fa — 7T, 0, 7r) and (0, — n, fa) (in Xi~ coor dhiates). The 
running coordinate \ = Xi or X — X3 is equal to fa (up 
to a trivial shift) and thus related to the external bias 
flux $ via x = — 27T$/$ - The current I = —cd^E then 
is given by the expression I.j(x) = {2e/H)d x Vf(x) and is 
double periodic in the interval < $ < $o ( we define 
the charge of the electron as — e and e > 0). Alterna- 
tively, one may measure the frequencies u>i^(x)i c f- Eq. 
directly via the resonant absorption of an ac-signal. 

At nonzero temperatures (but still T <C Ej) the in- 
duced potential is driven thermally and involves the free 



F i { X ,T)=Y J 



Tlog 1-e 



fi"i(x)/T 



(8) 



Thermal fluctuations become relevant for temperatures 
T > hjji and lead to an increase in the barrier 5Ff(T) = 
F t (*/2,T)-F f (p,T), 

8F { (T) = (fe*/2)[V2-(l + 2fioM)] 

+ Tlog[wf/2fi ]. (9) 

As a result, rather then decreasing, the phase stiff- 
ness in the tetrahedron increases with temperature and 
hence the Josephson current I,j(T) cx <9$i*f ($, T) in- 
creases with temperature until the fluctuation-induced 
potential disappears due to level broadening: ther- 
mally induced quasiparticles produce a level broaden- 
ing h/RC ~ (EjE c /hA)exp(-A/T), where we have 
used the Ambegaokar-Baratoff relation I C R ~ A/e. This 
broadening should remain small on the scale of the level 
separation Hu>{, from which we obtain the condition 
that T < A/ln(ftwf/A). Beyond this temperature, the 
Josephson current is expected to decrease again, result- 
ing in a non-monotonic behavior of Ij(T). 

In order to arrive at a fully quantum mechanical de- 
scription of the tetrahedron, we have to account for the 
tunneling processes between the points Oi, cf. Fig. 2. It 
is important to note that each set of 4 mid-edge points 
residing in one plane fa = has to be identified with 
one quantum mechanical state |Oi); on the other hand, 
the pair of classically degenerate lines connecting two 
states \Oi) and \Oj) in one of the faces describe differ- 
ent tunneling trajectories, which have to be added co- 
herently in order to arrive at the tunneling matrix ele- 
ment between the two states (the other two trajectories 
on the opposite face are equivalent). Let us concentrate 
on the tunneling process between \0\) and IO2); the two 
tunneling trajectories follow the lines parametrized by 
XUX2 € {0,-tt}, X3 € [0,7r], see Fig. 2(a). The tunnel- 
ing processes described by these two trajectories flip the 
phase across the four junctions '31', '01', '02', and '23' by 
±7r, which corresponds to a fluxon $0/2 traversing the 
tetrahedron as shown in Fig. 2(b) (the same arguments 
apply to the connected tetrahedron). The phase differ- 
ence between the two trajectories then corresponds to 
taking a full fluxon $0 around the two islands '1' and '2', 
which translates into the Aharonov-Bohm-Casher phase 
exp[27ri((7i + 32)]) with qi the charge on the island 'i' 
measured in units of 2e (see Ref. [22 for a detailed dis- 
cussion of charge-induced interference effects in small su- 
perconducting structures; these phases are generated by 
the topological term ft^i li^i i n the Lagrangian C in the 
presence of charges 2eqt, cf. the note below J3J)). Com- 
bining this phase factor with the modulus \a\ we arrive 
at the tunneling amplitude 



tu = -2|a| cos[7r(gi + q 2 )\ 



(10) 
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between the states \0\) and [Ox); a similar analysis pro- 
vides the amplitudes for all the other pairs. In the ab- 
sence of any charge frustration (i.e., for integer charge 
qi on each island) the system gains energy from hopping 
and hence t < 0, thus defining the sign in lfTU|) . The 
modulus |a| of the tunnelling amplitude follows from the 
semi-classical description of the one-dimensional motion 
under the barrier Vf(x) as given by Eq. and takes 
the form \a\ w [h/T(Q )] exp[-Sf (O )], with T(fi ) the 
classical period of motion and S[(Qq) the dimcnsionlcss 
action 2 ^, both evaluated at the ground state energy MIq, 



Si = 



1 



T = 



/32/ v v ■ r ■ 

V E C / Jo 



dx [\/2cos(x/2) - 1 - 2ft /^f] 1/2 



(Ej/Ec) 1/4 » 1, (11) 



Ec 

8h ( MV'" /J5j 



V 2 / „ x 2/3 



(12) 



where xo = 2arccos[(l+2f2o/wf)/v / 2]- The semi-classical 
analysis describing tunneling at the correct ground state 
energy Ml gives very accurate results, see below; the 
simple asymptotic form Ijlll) describes tunneling from 
the bottom of the well and is valid only deep in the 
quasi-classical regime. The result obtained here for the 
tetrahedron is smaller than the usual tunneling action 
S oc (Ej/Ec) 1 ^ 2 of a Josephson junction device involv- 
ing only the square root of the parameter Ej/Ec] the 
unconventional dependence on the ratio Ej / Ec is a con- 
sequence of the fluctuation-induced (rather than classi- 
cal) barrier and puts less stringent requirements (e.g., 
with respect to the smallness of the junction) on the fab- 
rication process of this new type of qubits. 

The above tunneling action can be probed via a mea- 
surement of the current-voltage characteristic of the de- 
vice: applying a fixed bias current across two islands, a 
finite voltage appears due to quantum and thermally in- 
duced phase slips. At low temperatures the appearance 
of quantum phase slips involves the finite action S ps = 
2hS{ and results in an exponentially small resistivity 
R tx exp(— 2Sf). At higher temperatures T > U ps h/S ps 
the phase slips are created thermally via activation over 
the barrier U ps = SFf. The phase-slip induced resistiv- 
ity then exhibits an unconventional saturation at high 
temperatures: increasing T beyond hwf{Ec / Ej) 1 / 4 , the 
resistivity first increases with temperature. As T > hwf 
the barrier SFf itself increases linearly in T and the fur- 
ther rise of R saturates at a value R cx Qo/cof; hence care 
has to be taken not to confuse the saturation in R at 
high temperatures with the (exponentially small) quan- 
tum resistance surviving as T — ► 0. 

Let us return to the quantum description of our tetra- 
hedron and study its low-energy spectrum as determined 
by the quantum coherent oscillations between semi- 
classical states \Oi); the appearance of the island charges 
2e qi in the tunneling amplitudes tij manifests itself in 
this level structure. Assuming a uniform distribution of 



the total charge 2e q on the isolated tetrahedron, the ma- 
trix elements take the form = t = —2\a\ cos{-Kq/2) and 
we have to diagonalize the matrix 



Ho = 




(13) 



Depending on the value of q, the low-energy spectrum of 
the isolated tetrahedron splits into singlets and doublets 
involving the energies E s — 2t and = —t, or remains 
tri-fold degenerate with E%(t = 0) = (here, energies are 
measured with respect to the unperturbed value — 2Ej + 
h(ujf/2 + f2o)); the ground state is a 



singlet if 
doublet if q 
triplet if q 



q = 4k, 
-4k + 2, 
■ 2k + 1, 



even, 
even, 
odd. 



(14) 



In the connected tetrahedron of Fig. 1(b) the charge 
is not quantized on the inner islands; however, the de- 
pendence (|10|l of the tunneling amplitude on the island 
charges qi remains valid and the above spectrum is recov- 
ered under appropriate biasing of the inner islands with 
either zero, one-quarter, or one-half Cooper-pair. Bias- 
ing the tetrahedron into the charge state qi — 1/2 then 
establishes a ground state doublet with eigenstates 



I-) 



= [|0 1 ) + C|Oa> + r|0 3 )]/V3, 

= [\o 1 ) + c\o 2 ) + c\o 3 )]/V3, 



(15) 



suitable for the implementation of a quantum bit; here, 
C = exp(27ri/3). These degenerate states (|15|) involve 
bonds resonating with opposite chirality in the device, 
see Fig. 2(a), and are reminiscent of the resonating dimer 
bond s in th e topologically protected qubits discussed in 
Refs. lirillL Also, such chiral states appear as degenerate 
spin-singlet ground states in the vanadium tetrahedron 
of the pyrochlore system^. The doublet |±) is protected 
by the gap A<j = 3t, separating it from the next excited 
(singlet) state 

|0) = [Id) + \0 2 ) + |0 8 )]/V3. (16) 

Combining the results i|l(J[l , (|ll|l , and l|12|) we obtain the 
protective and operational energy scale t of the qubit, 



t _ 1 / CI 
~E~i~4 \j2tt 



1/2 



Ec 
Ej 



2/3 



exp[-S f (n )]. (17) 



C. Charge limit 

We briefly extend our analysis to the charge-dominated 
regime with Ec 3> Ej. It is convenient to go over to 
a Hamiltonian description; starting from the above La- 
grangian (cf. Eqs. and ©) and eliminating the vari- 
able 4>q one obtains the expression 



H 



EclQl 



Qi 



Ql + (Qi + Q2 + Q3 



(18) 
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where Qi — ~id ( j >i — qi and qi is the induced charge 
on the i-th island (in units of 2e), qi = J^jCijVj for 
the connected tetrahedron, while an additional term 
[q - (C/A)J2j v o\l^ with C = Ey C H the total capac- 
itance, has to be added for the isolated tetrahedron 22 
(here, dj denotes the capacitance matrix, see ©, and 
Vi are the bias potentials applied to the islands) . For the 
isolated tetrahedron the total charge q — J2o Qi i s integer, 
while the total induced charge q = J^i Qi can take any 
value in the connected device. The Hamiltonian 118(1 de- 
scribing both devices becomes identical under symmetric 
bias and for specific values qi = k/4 with k an integer; 
under these conditions the maximal symmetry S4, i.e., 
the tetrahedral symmetry Td, is established. Note that 
a symmetric bias with equal charges qi = g/3 on the 
three inner islands of the connected tetrahedron in gen- 
eral guarantees only for a S3 symmetry. 

We determine the spectrum for the uniformly charged 
isolated tetrahedron with q = 4k + 2. In the limit 
Ec S> Ej the operators Qi take on integer values and the 
charging term is minimized by distributing two bosons 
onto the four sites avoiding double occupancy (in this 
limit, the term E c (Qi + Q2 + Q3) 2 in ifTFIl describes the 
charging energy EcQ\ of the middle island). The result- 
ing six states '01', '02', '03', '23', '31', and '12' are degen- 
erate with an energy Eq = 2Ec- The hopping term 
lifts this degeneracy through the mixing via the Joseph- 
son coupling Ej: each state 'zj' hosting Cooper-pairs on 
the islands H' and 'j' exchanges particles with all other 
states except for l kV, where k, I ^ i,j. The Hamiltonian 
describing the mixing of the six two-Boson states may be 
written as a matrix product 



H 2B = It t J ® ^ 1 J, (19) 

where the sign of the tunneling amplitude t = Ej /2 > 
is positive for our frustrated tetrahedron (again, we 
choose a symmetric gauge with all Josephson couplings 
reversed in sign). The eigenvalues of the direct matrix 
product in l|19l) are given by the product of the eigenval- 
ues (2t, —t, — i) of its first factor and those of the second 
factor, (0, 2); correspondingly we find the first 6 levels at 
the energies 

E d = -Ej doublet, 

E t = triplet, (20) 

E s = 2Ej singlet. 

The first excitation now is a triplet rather than the singlet 
state found in the opposite limit Ej 3> Ec', hence, de- 
creasing Ej/Ec from large values, the singlet and triplet 
energies cross each other and the first excitation gap 
changes over from 3t (at Ej ^> Ec) to 2t = Ej (for 
Ec ^> Ej). The precise location where this crossing ap- 
pears can be found from the numerical analysis described 



below. The same analysis can be repeated for the con- 
nected tetrahedron; the six states degenerate under the 
capacitive term then involve either one or two charges 
on the three inner islands and the mixing term V v de- 
scribes charges hopping between the three islands as well 
as hopping of one charge to and from the ring. 



D. Numerical Results 

The above results can be verified numerically via di- 
agonalization of the Hamiltonian l|18l) with the help of a 
Lanczos algorithm; in the charge basis the mixing term 
Vir then describes the hopping of charges between the is- 
lands. Going to a phase representation, the bias charges 
qi are conveniently accounted for via the boundary condi- 
tion for the wave function, \f r n ((/>i+27n5ifc, 02 + 27T(52fc, </>3 + 
27r(5 3fe ) = exp(-27rzgfe)\P„(</>i,(/>2,</>3), k = 1,2,3, after 
a suitable gauge transformation 2 ^. The results of such 
an analysis for the charge state qi — 1/2 is shown in 
Fig. 3, where the excitation gap Ad protecting the qubit 
states against higher excitations is shown as a function 
of Ej/Ec- The crossover from the charge to the phase 
dominated regime, where the singlet and triplet excited 
states cross, takes place at Ej/Ec ~ 5. The analytic re- 
sults (|17fl and (|20|l describe well the data away from the 
crossover regime. One expects the quasi-classic result to 
become exact in the limit of large Ej/Ec; however, the 
result lfT7|l has been calculated using the one-dimensional 
approximation Vf for the potential and one cannot ex- 
pect perfect agreement with the numerical data. Still, 
the quasi-classic approximation turns out accurate over a 
very large regime extending down to parameters Ej / Ec 
of order 10: scaling the dashed line in Fig. 3 by 0.8 the 
quasi-classic result cannot be distinguished from the nu- 
merical data for Ej/Ec > 20. 

The numerical results show that the suppression of the 
tunneling amplitude t is indeed small in the tetrahedron. 
E.g., choosing a value Ej/Ec of order 100 the energy 
scale t of the qubit is suppressed by a factor ~ 2/1000 
with respect to the energy scale Ej of the junctions. For 
a conventional device this suppression involves an action 
S « C\J Ej/Ec, with the numerical c depending on the 
specific setup. E.g., for the 4-junction loops studied in 
Ref. |3 the numerical era 1.6 and choosing the same 
value of Ej/Ec this implies a suppression of quantum 
fluctuations by a factor ~ exp(— 16) ~ 10~ 7 ; such a de- 
vice then resides deep in the semi-classical regime and 
quantum effects are heavily suppressed. 

In summary, we can tune our tetrahedral structure 
such as to realize a doubly-degenerate groundstate corre- 
sponding to a spin-1/2 system in zero magnetic field; the 
device can be realized using moderately large junctions 
with Ej/Ec of order 100 while keeping an appreciable 
operational energy scale t, a consequence of the partic- 
ular frustration in the device. Below, we will show that 
this ground state remains robust to quadratic order in 
the external noise and hence provides a suitable starting 
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FIG. 3: Excitation gap Ad protecting the qubit doublet 
against higher excitations as a function of Ej /Ec; the dashed 
line is the semi-classic result based on 1)17^ . The inset shows 
an expanded view illustrating the crossing of the singlet and 
triplet excited states as the system enters the charge dom- 
inated regime at low values of Ej/Ec; the dashed line is 
the analytic result I20H valid in the charge dominated regime 
Ec S> Ej. Data obtained with help of a Lanczos algorithm 
operating in the charge representation with 27 charge states 
between q = ±13. 



point for the construction of a qubit; we will concentrate 
on this doubly-degenerate case and its use for quantum 
computing in the following. 



E. External fields 

Fabrication errors and external bias induce splittings 
and shifts in the levels. With respect to the qubit's func- 
tionality, random external signals produce decoherence, 
while prescribed bias signals are used for its manipula- 
tion. Fabrication errors mainly affect the coupling Ej of 
the junctions — we denote the corresponding deviations 
by eijEj (note that the fabrication tolerance is improved 
for larger junctions, leading to smaller values of ■). Ex- 
ternal bias signals appear randomly through fluctuations 
in the local magnetic field and through stray charges; ar- 
tificially generated signals can be applied through prop- 
erly placed current loops biasing the sub-loops of the 
tetrahedron (currents ii in Fig. ^) and through capac- 
itive charging of the islands (via voltages Vi). We denote 
the corresponding bias fields by Sf = — 27r$i/$o an d 
Sf = Sqi, i = 1,2,3. E.g., for the isolated tetrahedron 
the flux $3 penetrating the sub-loop '0-1-2' is encoded 
in a bias angle <5f = S((f>2 — <j>\) = — 2-7T $ 3 /<E>o along 
the link '1-2'; equivalent definitions apply to the other 
sub-loops with the cyclic replacement 3 — » 1 — > 2 — > 3. 
For the connected tetrahedron we have to account for 
the total flux <E> threading the outer ring; assuming a 
symmetric setup, this flux will induce the phase shifts 
8 = — 27r$/3"I > o m each of the segments 'ml-m2', 'm2- 
m3', and 'm3-ml', cf. Fig. 1(b). The bias angles Sf along 



the links l j-k' then have to account for this flux via the 
modified form Sf = 2tt ($ 4 - $/3)/$ , where $i denotes 
the flux through the loop 'mj-mfc-fc-j-mj'. Note that we 
have assumed that the currents circulating in the outer 
ring are still sufficiently small such as not to produce sig- 
nificant self-fields. While the charge bias Sf only affects 
the tunneling matrix elements £y, i.e., the kinetic energy 
term, the flux bias Sf modifies both the potential and 
the kinetic energy terms in the Lagrangian. 

We first determine the modification of the potential 
energy V n , cf. JIJ, due to applied fluxes Sf, 



2E, 



cos [xi 



x cos xi 



(- - f ) 



(» + f)cos(»-f)+co S ( X3 + f) 



COS X2 



5J 
2 



COS X3 



SJ 
2 



(21) 



Assuming small perturbations, we expand (|21l) in Sf <C 
Uf/Ej] we concentrate on the point O2 — (0, — 7r, 0) and 
combine the result with the kinetic term to arrive at the 
Hamiltonian JfJJ with the additional term 



SH 
~E~j 



$i($X3 - 5x2) + <5f (<$X3 - Sxi) 



+ St(S X 2-S X i) 



\{Sf + 5f 



(22) 



Classically, the force term in l|22|l lowers the energy in- 
definitely as the system runs away along the degenerate 
classical minimal lines; e.g., a perturbation <5| > pro- 
duces a runaway either along the (0, 0, 8x3) direction or 
along the (—8x1, 0, 0) direction. However, quantum fluc- 
tuations generate a finite potential along these lines, re- 
sulting in a linear response in the coordinates Sxi and 
a quadratic change in energy v%. Indeed, second-order 
perturbation theory in the force term of (|21|l produces 
the result 



V2_ 

Ej 



(Sf 



($2 



sir 



+ (eo2 — £03 — e 0l) + ( e 31 — e 12 ^ £23)1 (23) 



with 



1.0 (Ej/Ec) 



1/3 



(24) 



obtained from a numerical solution of the perturbed 
eigenvalue problem combining © and l|22|l . In (|23|l we 
have dropped the term (Ej/2)(Sf 2 + <5f 2 - 5| 2 ) as it 
is small by a factor (Ec / Ej) 1 ^ 3 compared to the lead- 
ing term; also, we have added a term due to deviations 
tijEj in the junction couplings. Equivalent expressions 
for the other minima follow from cyclic permutation of 
the indices. As a result, we obtain the relative shifts 

V ij = V j ~ V i (h3)k = MA k ^ t,j) 



2(eoj — eoi 



3 



251 



(25) 
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in the minima. First, we note that small random fluxes 
do not affect these positions in the first order of these 
fluctuations; the corrections appearing in quadratic or- 
der then are small. Second, we note that fabrication 
errors in the junction couplings can be compensated 
by appropriate choices of bias fluxes Sf . 

In the determination of the perturbed tunnneling ma- 
trix elements Uj , we ignore the modifications arising due 
to fabrication errors and concentrate on the effects of 
flux and charge signals, random or externally applied. 
By way of example, we calculate the tunneling matrix 
elements ti2 and t2i connecting the states \0±) and |Oa)- 
The presence of perturbing fluxes Sf and 5 2 shifts the 
potential V n by v along the line 7, to lowest order in 
Sf, v(x) = -Ej sinx(5* + £f)- This shift produces the 
changes SSf 2 = Ts(Sf + S 2 ) (the correction <5S^ applies 
to the trajectory 7 in Fig. 2(a)) with 



this calculation into the form 

-Hqubit = e d 1 + h • o 



(30) 



with o = [a x , a y , a z ) the Pauli matrices; switching of the 
'magnetic' fields h x and h z then produces the standard 
amplitude- and phase-shift operations required for the 
manipulation of the individual qubit. The shift ea and 
the effective 'magnetic' field h are given by the expres- 
sions 



s w 1.5 (Ej/Ec 



\3/4 



(26) 



ed 
h x 

hy 

h z 



-[vi+V2 + v 3 -(tt 2 + t s 23 +t s 31 )}, 

- [vi -~{V2+V 3 ) + (2^ - t\ 2 - t\ x ) 

^-^) + ^(*5 a -«x). 

— ~/=(*12 + *23 + *3l)i 



(31) 



in the action Si determining the modulus |a| of the tun- 
neling amplitude (JTUJ; as before, the expression lj2~6l is 
valid deep in the semi-classical regime. 

The presence of perturbing charges Sf and S 2 modifies 
the Aharanov-Bohm-Casher phase associated with the 
two trajectories (a charge Q encircling a flux $ counter 
clockwise produces the phase exp[27ri($/$o)(Q/2e)]): 
they pick up the additional phases exp[±iir(8f + 6®)}, 
with the plus sign belonging to the trajectory 7, cf. 
Fig. 2(b). Combining the perturbations in the fluxes 
and charges, the change in the tunneling amplitudes 
5ti2 — t\ 2 + it\ 2 = St 21 takes the form 



Stu 
t 



; (27) 

expanding the exponentials, the symmetric and antisym- 
metric parts are given by the expressions 



t{ 2 = t[ S 2 (Sl+8*) 2 /2-Tr 2 (8?+5f) 2 /2], 



(28) 



further terms quadratic in <5* (e.g., arising from the next 
term in the expansion of Sf) are small by the factor 
(Ec/Ej) 1 / 4 . With respect to the qubit's stability, we 
note that all perturbations appear in second order of the 
small quantities Sf and Sf . In the further analysis below 
we will drop the flux bias term cx Sf 2 in tfj against the 
(parametrically large) shifts m in the potential energy. 

We determine the new energy levels perturbatively: 
the perturbation SH, written in the space of semi- 
classical ground states |0,), takes the form 



vi 



1-12 



a 12 i 3 i 



SH 



L 12 



n 12 
<r 31 I 23 



»2 



*23 + 



iE 31 
^23 



(29) 



^23 



l<3 



Next, we find the corresponding matrix elements 
(±\SH\±) in the projected space spanned by the dou- 
blet |±), cf. i|15|) . It is convenient to cast the result of 



note that the perturbations Vi in the potential come with 
the large amplitude Ej, while those in the kinetic energy 
(t s ' a ) involve the smaller energy scale t of the tunnel- 
ing matrix element; in (|31|l we keep both terms as we 
might choose to manipulate the qubit via changes in the 
charges S Q alone. Expressing the perturbations in terms 
of the flux- and charge-parameters Sf and Sf we obtain 
the results 



ed = 



-^p (sf + sf + sf + <s?aj + sfst + 

t-K 2 

m + si? + & + s*) 2 + + s?n 

vE 



j 



hy = 



6 

tix 2 

+ T 

vEj 
2V3 
tir 2 



m+siy-(st+s^y-(sf+si) 



[2Sf - Sf - Sf - 2S§{5% ~ Sf) 



- 2^ 



S?)}, 



(28f 



2^3 

7T St 

7s 



[Sf-Sf + 2S?(S$-S$)}, 



(32) 



+(S% + 6Z)(6° + 5$) + (SI + 5!)(5* + S?)}. 



For convenience, we also cite the energy shift e s of the 
singlet state, 



1 



[vi+v 2 + v 3 + 2{t s 12 +t s 23 +t s 31 )], 



(33) 



3 v 1 



Sf + St5* 2 +8t5t + 8*5t) 



x m + s c i? + (5f + s«) 2 + + s?n 

While the bias variables Sf and Sf as well as the fields ed, s 
and h z do respect the symmetry of the tetrahedron, the 
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bias fields h x and h y do not. The tetrahedral symmetry 
may be made explicit by reexpressing the planar field hj_ 
in the hexagonal basis ei = (1,0), e 2 = (— 1/2, — %/3/2), 
e 3 = (-1/2, v/3/2), h ± = J2i Mi with 

hi = ^ (SI + 5tf+ X 2«?(J? + 5«) 2 ]; (34) 

the relation between the tetrahedral symmetries and the 
rotations in spin space then becomes obvious. The three 
symmetric fields hi relate to the two cartesian fields via 
h x = h\ — (/12 + hs)/2 and h y — V3(h 2 — h 3 )/2; there is 
no inverse transformation. 



components h x , y — h±exp(iu>h ± t), while keeping the 
z-component h z fixed, defines the operator J/ecrry = 
expier z f2/2 after one period of rotation. Here, f2 = 
27r(l — h z / \Jh 2 z + JvjJ is the solid angle spanned by the 
rotating field-cone. Selecting a field vector h pointing 
at an angle of 60° with respect to the z-axis, the result- 
ing operator exjp(ina z /2) changes the relative sign of the 
components along |±), i.e., the eigenstates \s,a) of the 
operator a x transform into one another upon each pe- 
riod of field rotation. Note that operator E/serry does not 
depend upon the rotating field frequency as long as 
w^ 1 is much shorter than the qubit's decoherence time 

^-dcc ■ 



III. NOISE SENSITIVITY AND 
MANIPULATION 



IV. MEASUREMENTS 



The above results exhibit two remarkable features of 
our tetrahedral qubit: 

i) Both, charge and flux noise appear only to quadratic 
order, a feature which is easily traced back to the absence 
of polarization charges and currents in the qubit's ground 
state. This guarantees for a long decoherence time of the 
tetrahedral qubit, similar to the 'quantronium' discussed 
by Vion et aZi Numerical analysis confirms the weak 
susceptibility to charge noise: a uniform random bias V" € 
[—Vq, Vq] acting on the islands produces a small doublet 
splitting S » 0.2 Ad(2eVb/^c) 2 , roughly independent of 
Ej/Eq- At the same time, the quadratic dependence 
on bias allows for a qubit manipulation via ac-fields and 
the dangerous low-frequency noise can be blocked with 
appropriate filters. 

it ) The tetrahedral qubit admits a large variety of ma- 
nipulation schemes using either magnetic or electric bias. 
The 'planar fields' h x and h y can be manipulated via 
changes in flux alone, e.g., setting S 2 = <5f = 5* ^ Sf/a x 
with a x = (— 1 ± 2), we can direct the field along 



the x-axis, while the flux state 8f — 8*1(1 + a 2 )/2] 



1/2 



S* , with a y = — (2 + Vo) produces a field 



pointing in the y-direction (here, we assume no charge 
bias, Sf = 0; with these choices of parameters the sub- 



leading term Ej(8f 



£f )/2 in (23 contribut- 



ing to h x with (Ej/6)[8$ 2 + <5f 2 - 25 f] and to h y with 
{Ej/2V3)[8f - 5f], vanish as well). Proper choice of 
amplitudes and phases allows for the generation of a ro- 
tating planar field. The axial field h z , however, involves 
a modification of fluxes and charges, cf. (|32|) ; choosing a 
uniform flux- 5f — 5* and a uniform charge-bias Sf = S Q 
produces a pure axial field 



h, = -4irV3stS' s '5 c ' 



(35) 



The free manipulation of the 'magnetic field' h(t) 
allows for the implementation of Berry-phase type 
phenomena^SiSL 2 ^. E.g., the following provides a sim- 
ple realization of the NOT-operator in the basis \s,a) = 
[|+) ± |— )]/y2: adiabatically rotating the transverse 



Finally, we discuss several potential procedures for the 
measurement of our qubit's state. We assume an ideal 
symmetric device. As the qubit reacts to external bias 
only in second order, the measurement of its state in- 
volves a two-step process: In a first step, the qubit 
is pushed away from the symmetric point of operation 
through appropriate driving via external charge (Sf) and 
flux (8f) bias. As a consequence, the previously quiet 
qubit state now develops finite internal currents and po- 
larization charges. In a second step, these signals have 
to be measured by appropriate devices and the meter 
reading will tell about the qubit's state. 

As a simple illustration we consider the measurement 
of the qubit's state in the basis |±), i.e., the operator a z , 
using a charge bias 8f <C 1. This charge bias generates a 
current flow within the qubit and the associated flux can 
be measured by an external SQUID loopS, cf. Fig. Ufa). 
We express the matrix (|2*9l in the basis {|+), |— ), |0)}, 



SHo = 



ed 



ih y e 



ih, 



Oy 



h x + ih y e d - h z 



e ih Qy 



(36) 
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ih, 



0y 



with ed, s; h x ,h y ,h z given in (|31|l and (|33|) above and 

*3l)] 



4 = 7 [2v! - v 2 - v 3 - (2t s 23 - t s 12 



± 



2\/3 



(2^23 t 31 ), 



(37) 



h± - 
0y 2^3 



A finite 'magnetic field' h = (0, 0, h z ) induces the shifts 
8E ± = ±h z in the two states |±). Charge-biasing the 
device induces the current (we assume Sf -C 1) 



2e dSE^ 



= T- 



n dSf 

2-Kest 



(38) 



Vsh 
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in the corresponding sub- loop '0 — j — fc' (the above signs 
apply to the isolated tetrahedron and have to be reversed 
for the connected device; a positive current runs counter 
clockwise around the loop). Alternatively, applying a 
(small) external flux-bias induces the voltage (see Fig. 
Btb)) 

V? - ~ (39) 



2e 88? 

TTSt 



St) 



on the island H 
tron transistor 



which can be measured via a single elec- 
(SET) device^; here i,j,k G {1,2,3} 
and pairwise different. Choosing the bias in the ratio 
g*,Q . g*,o . s *,q = 3 : i:i H m its the current /voltage 

to the loop/island H'. 
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FIG. 4: Measurement setups allowing the identification of the 
qubit state, (a) An external charge bias induces currents in 
the qubit structure with a circularity depending on the qubit 
state. The flux associated with these currents is inductively 
coupled to a SQUID which is driven close to criticality during 
the measurement. Depending on the qubit flux, the SQUID 
is driven overcritical and the associated voltage is measured, 
(b) An external current bias induces polarization charges in 
the qubit structure with a polarity depending on the qubit 
state. The associated charge is capacitively coupled to a SET 
which is driven close to (charge) frustration during the mea- 
surement. Depending on the qubit's polarization charge the 
SET is driven into the conducting state and the associated 
current is measured, (c) Similar to (a), but with qubit cur- 
rents directly channelled through the measurement junctions 
with large couplings _B m Ej; again, the presence of a volt- 
age Vm in the external loop carries the information on the 
qubit state, (d) States of the measurement junctions (with 
phase 7) at fixed classical driving current 7 m with a slowly 
decaying state |0) (decay rate To) and a fast decaying state 
|1) (decay rate Fi S> To) depending on the qubit state. We 
assume a slow energy relaxation for the qubit, Te <C Y\. 

Several issues have to be considered in the measure- 
ment process, cf. the discussion in Ref. l3ll Relevant pa- 



rameters are the dephasing- and mixing rates T v and T e 
induced by the measurement apparatus and their relation 
to the separation Eqi between the qubit eigenstates; in 
a weak measurement scheme we have T v <C Eqi, while a 
projective measurement is characterized by a strong cou- 
pling with Eqi <C T v such that the quantum evolution of 
the system is quenched rapidly. While dephasing trans- 
forms a coherent superposition of states into a classical 
mixture through elimination of the off-diagonal elements 
in the qubit's density matrix, the mixing induces transi- 
tions between the qubit's states and thus spoils the mea- 
surement. Usually, a good measurement setup makes use 
of decoherence in transferring classical information to the 
measurement device but avoids mixing, hence T v 3> Te- 
The simplest situation is realized when the measured 
observable commutes with the qubit Hamiltonian — in 
this case the measurement preserves the qubit's eigen- 
states. On the other hand, if the measured observable 
does not commute with the qubit Hamiltonian, the mea- 
surement has to be completed before mixing sets in; 
hence < r mcas < T v <C Eqi/H for a weak measure- 
ment, while the sequence Te < r moas < Eqi/H <C T v 
applies to the projective measurement (here, r mcas = 
1/^meas denotes the inverse measuring time). 



(b) 






FIG. 5: Current-voltage characteristics of meter devices. The 
solid/dashed lines refer to the different quantum states of 
the qubit shifting the characteristic of the meter device, (a) 
In an overdamped SQUID the current-voltage characteristic 
is single valued, with a voltage due to a finite rate of in- 
dividual phase-slips. J m is the imposed measuring curent. 
(b) The current in the SET results from a continuous flow 
of individual electrons traversing the island. V m is the im- 
posed gate voltage during measurement. (c) An under- 
damped SQUID/ Josephson junction exploits an instability 
where a single phase-slip triggers a transition to the dissi- 
pative branch. 

The measurement schemes described in Ref. |^ involve 
dissipative meter devices, e.g., an overdamped SQUID 
with a well defined (single-valued) current- voltage char- 
acteristic as depicted in Fig.[S{a) or a SET with the char- 
acteristic shown in Fig. EJb) . In both cases the measure- 
ment involves numerous dissipative events, either many 
phase-slips producing the voltage across the SQUID de- 
vice or many electrons traversing the island of the SET. 
The fluctuations due to the phase slips/electrons act back 
on the qubit, enforcing its loss of phase coherence. Such 
measurement schemes can be implemented in terms of 
weak or strong (projective) measurements. 

This type of measurement has to be contrasted with 
a meter characterized by an instability, such as an un- 
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derdamped SQUID or Josephson junction with a char- 
acteristic as shown in Fig. OJc). Such a meter does not 
couple dissipatively to the qubit but switches to a dissi- 
pative state only after the measurement, e.g., after the 
occurrence of one phase slip, hence T v <C Eqi; the mea- 
surement is weak and generates a small imaginary part 
to the energy of the qubit eigenstate which then may be 
determined in a decay process. The qubit states are iden- 
tified through their (exponentially) different decay rates. 
Such a measurement scheme has been used recently by 
Vion et ali. In their setup, an additional measurement 
junction is introduced into the qubit loop. The current 
generated by the qubit is superimposed on an external 
measurement current and drives the measurement junc- 
tion towards criticality, see Fig. Hfc). Here, we make 
use of the symmetric setup shown in Fig. E{a) involv- 
ing six classical measurement junctions with equal cou- 
plings E m 3> Ej . The measuring currents J m j are fed into 
the device through the points 'nl', 'n2', and 'n3' and re- 
moved at the points 'ml', 'm2', and 'm3'. The measuring 
current I m flowing in the segments 'n3-ml' and 'n3-m2' 
(and in the other obtained through cyclic permutation) 
is chosen close to the critical current J m c = 2eE m /H of 
the measurement junctions, requiring to feed a current 
7 ext » 2I m>c into the tetrahedron via the input lines 'nl', 
'n2', 'n3' and retraining them via the output lines 'ml', 
'm2', 'm3'. A convenient measuring setup is the symmet- 
ric one with equal currents crossing the junctions ! n3-ml' 
and 'n3-m2' (and other pairs obtained through cyclic per- 
mutation). Driving the qubit with external bias fields 
8 1 and 8f induces additional currents I^i in the loops 
'mj-mfc-fc-j-mj' which are characteristic for the quan- 
tum state of the qubit. Depending on the relative flow 
direction of the qubit currents and the measuring cur- 
rent, the six measuring junctions are either driven to- 
wards (I m + I Chi ~ J m)C ) or away (J m + I q>i < I m , c ) from 
criticality. The voltage pattern appearing on the six junc- 
tions then allows for the identification of the qubit state. 
Within this scheme, care has to be taken not to spoil the 
symmetry of the measurement setup by a circular cur- 
rent flowing in the outer ring, hence the system should 
be flux biased such that the total flux $ through the ring 
remains zero. This can be easily achieved by compensat- 
ing the fluxes $i in the loops 'mj-mfc-fc-j-rnj' through 
an equal and opposite flux $a = — 2<=i ^* through the 
central triangular loop. 

Below, we study two schemes, a weak measurement in 
the Hamiltonian basis with T v < Eqi and a projective 
measurement onto the current basis with Eqi < T v . In 
both cases we discuss the measurement of two operators, 
the measurement of a z and of a x . 



A. Measurements in the qubit's eigenbasis 

A proper measurement in the energy eigenbasis im- 
poses a number of constraints on the measuring device. 
First of all, the phase decoherence rate has to be small, 




FIG. 6: (a) Symmetric setup measuring the quantum state of 
the qubit. Currents below twice the critical current / m . c = 
2eE ul /h of the six equal measurement junctions are fed into 
the qubit loop through the points 'nl', 'n2', and 'n3', and 
removed through the leads at 'ml', 'm2', and 'm3'. The 
qubit current / q flowing in the segments 'n3-ml' and 'n3- 
m2' (and cyclically permutated) add to these measurement 
currents and drive the junction towards or away from crit- 
icality. The voltage pattern measured on the six junctions 
tells about the qubit's state. Care has to be taken not to 
spoil the symmetry in this measurement setup through a ring 
current in the outer loop; such a current can be avoided by 
appropriate compensation of the flux bias through the indi- 
vidual loops l mj-mk-k-j-mj' with an equal and opposite flux 
through the central loop '1-2-3-1'. This can be realized, e.g., 
via an independent control on the individual loops and on the 
overall flux penetrating the entire structure, (b) circuitary in- 
troducing damping on the six measuring junctions. Driving 
the parallel junctions (with couplings E'j) overcritical via a 
large current 27 switches the parallel channel into the resis- 
tive state. The contacts 'dl' and 'd2' of such a dissipative 
element are then connected to the points 'ml', 'm2', and 'm3' 
of the qubit structure in order to make the associated currents 
classical. 



T v <C Eq\. The decay rate T\ via tunneling of the high 
energy state (low barrier state, cf. Fig. 0Jd)) should be 
large, Ti 3> T^; such that the system tunnels before de- 
caying to the (metastable) ground state. Finally, the 
measurement time i mca s should be larger than the inverse 
tunneling rate, i mC as > . Note, that phase decoher- 
ence before tunneling, i.e., T v > Ti, is not a necessary 
requirement in this type of measurement; the decay may 
as well proceed out of the coherent state, i.e., a superpo- 
sition of the qubit states. In this case, the first phase slip 
triggers the projection, while the subsequent phase-slips 
produce the large voltage signal. 

In order to measure the operator cr n , i.e., the spin pro- 
jection onto the axis n, we first apply a 'magnetic field' 
h = hn directed along n. This is achieved via proper 
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charge- and flux-biasing as described by the equations 
(J22J. The doublet space {|+), |— )} then is split with new 
qubit eigenstates |+) n and |— ) n separated by the 'Zee- 
man' energy AE = 2h. The currents 1^-^ flowing in 
the segments 'mj-nfc' across the measuring junctions are 
equal to the loop currents in 'mj-mfc-fc-j-mj' and 
thus are determined by the derivatives 



In 



2e dh 
JidS? 



(40) 



In order to avoid the flow of a circular current in the 
outer ring, the driving flux bias 5f has to be properly 
compensated as described above. 

We proceed with the evaluation of the currents H4U|) 
associated with the measurement of a z and a y and gen- 
erated by the application of magnetic fields h z and h x 
along the z- and x-axis, respectively. Starting with the 
measurement of <r z , the projection along the z-axis, we 
choose a uniform charge and flux bias (— » h = (0, 0, h z ), 
cf. J22J)) and obtain the qubit loop currents 



q ' 1 V3» 



(41) 



i = 1,2,3, in the qubit state |±) z . Hence all loops are 
equally driven and voltage signals on the triple 'nl-m3', 
l n2-ml', l n3-m2' identify the qubit state |+), while finite 
voltages on the junctions 'nl-m2', 'n2-m3', l n3-ml' are 
associated with the |— ) state, provided that 8 Q > 0. 

Second, the projection along the x-axis can be mea- 
sured by applying a 'field' along h x using a flux bias 
$2 = S 3 = = Tnis bias generates a field 

h x = 4^_E/((P) 2 /3 and imprints qubit currents of mag- 
nitude 



7 , Sf^* 

q < 4 + 3h 



(42) 



for i = 2, 3 in the qubit state |±) x , while I* l =0. Hence 
the pairs 'n2-ml', l n3-m2', and l n2-m3', 'n3-ml' are equal 
and oppositely driven, while the junctions 'nl-m2' and 
'nl-m3' do not experience any additional drive due to 
the qubit structure. 



B. Projective measurement in the current basis 

A projective measurement in the current basis requires 
to switch on a strong decoherence T v ^> Eqi during the 
measurement. This decoherence then projects the state 
of the qubit at onset of dissipation onto the current ba- 
sis and keeps it there via the Zeno (watchdog) effect, 
cf. Ref. OOJ A suitable circuit allowing to turn on de- 
coherence is shown in Fig. EJb). The admittance be- 
tween the points 'dl' and l d2' is given by the expression 
Z~ x (w) = 1/2R - E' m /2ihu}R Q , with R Q = h/4e 2 the 
quantum resistance. Choosing parameters R < Rq (this 
guarantees a sufficient decoherence in the 'on' state, see 
below), E' m ~ (10 — 100)i?j, and an operating frequency 



w < O.lEj, we find that Z w -2iR Q (fvjj/E' m ). Hence, 
at zero applied current I the conductance between the 
points 'dl' and 'd2' is dominated by the Josephson junc- 
tions and is mostly imaginary at low frequencies, allow- 
ing us to ignore dissipation ('off' state). However, when 
the system is biased by a current I larger than criti- 
cal, the Josephson current disappears and the resistance 
R provides a significant source of dissipation as quan- 
tified through the dimensionless parameter a — Rq/R. 
In addition, while the Josephson junction itself involves 
a large quasi-particle resistance at low temperatures, 
when driven with a large current the junction switches 
to the resistive state with the resistance R' m ; the latter 
is related to E' m via the Ambegaokar-Baratoff relation 
R'^/Rq — A/E' m . The projective measurement in the 
current basis then starts with switching on a strong dis- 
sipation (such that 3> Eqi) which is diagonal in the 
eigenbasis of the qubit 's current operator J q . 

Before continuing with the discussion of an appropriate 
projective measurement in the current basis, we discuss 
one more issue related to the quality of such a measure- 
ment. The above measurement schemes (via coupling 
to a SQUID or a SET) are diagonal in the qubit's sub- 
space spanned by |+) and |— ) and can be implemented 
in terms of weak and strong (projective) measurements, 
cf. Ref. 1^1 However, extending the analysis to the low- 
energy subspace {|+), |— ), |0)}, we find for the current 
operator I\ = (2e/h)dSHo/dSf \s* =0 the expression (cf. 



2iTest 



i*?2 + %l) 
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where Sfj = Sf +5®; choosing a charge-bias Sf /3 = S Q 
—5® = — 5® this reduces to 




(43) 



Similar expressions apply to the other current operators 
1% and ^3 and, replacing 2Trest/Hy/3 — * nst/2eV3 in the 



prefactor and substituting Sfj — > 5fj = Sf 



Sf, to the 



voltage operators % = (l/2e)dSHo/dSf . 

Hence a projective measurement onto the current ba- 
sis is in fact non-ideal as the measured observable, the 
current, does not commute with the Hamiltonian when 
going beyond the qubit sector. In a 'high quality' mea- 
surement of our qubit state one would request that the 
off-diagonal matrix elements (between the doublet and 
singlet states) of the operator are much smaller than the 
matrix elements in the doublet subspace. Otherwise, the 
measurement has to be repeated many times in order 
to arrive at a proper readout of the qubit state. We 
then need to identify special measurement configurations 
where the off-diagonal matrix elements remain small or 
even vanish altogether. 
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We will now give a specific example how this goal can 
be achieved for a projective measurement of the rx z op- 
erator. We will show that a charge-bias with 



6? = 



(44) 



results in a high-quality measurement of a z and a y if 
we choose a large bias 5 Q = 1/2. This bias induces a 
current flow of equal magnitude but different chirality 
in the loops (cf. Fig. Eta)) 'ml-l-3-m3-ml' (current I2) 
and 'm2-2-l-ml-m2' (current I3; loop currents circulating 
counter clockwise are positive); no current is induced in 
the loop 'm3-3-2-m2-m3'. The total current through the 
link 'ml-1' is given by / (ml_1) = I 2 — I 3 and will be used 
to drive the measurement current I m fed symmetrically 
into the tetrahedron at the points l m' ovcrcritial. We first 
describe the setup projecting the link current p all ~ 1 '> and 
subsequently derive an explicit expression for the current 
operator. In a third step, we then describe the actual 
measurement process in detail. 

We will show below that the link current / (ml_1) will 
not mix to the third state |0). However, this property is 
not shared by the individual loop currents I2 and I3 . It is 
then important to devise a setup projecting the state onto 
the link current /( ml_1 > rather than the loop currents I2 
or ^3 . This is achieved by the symmetric coupling of three 
'dissipators' in between the points l ml-m2', 'm2-m3', and 
'm3-ml'. Turning on one of these dissipators renders the 
current in the corresponding link classical. In our setup 
we wish both loop current, I2 and ^3, to be classical and 
hence switch on the dissipators 'ml-m2' and 'ml-m3'. 
In this situation it is the link current /< ml_1 > which is 
coupled to the dissipative reservoir (the segment 'm2- 
m3' remains quantum and can be contracted to a point, 
cf. Fig.d) and hence the quantum state is projected to a 
state with fixed current J'" 1 - 1 '. In an ideally symmetric 
setup this current will be symmetrically split into the two 
arms 'ml-m2' and 'ml-m3' and, depending on the qubit 
state, will drive a specific set of junctions overcritical. 




dissipative ^ 
environment 



FIG. 7: Dissipative measurement setup projecting the qubit 
state onto the fixed link current /I™ 1 - 1 '. 

Let us then turn to the calculation of the link current 
/(mi-i) j n a gjgt s tep ; we have to generalize the tunneling 



matrix elements <|28[1 to allow for a large charge bias, 

t + % = tcosinSfj), 

t% = s^-sm^g)- ( 45 ) 

Assuming the above form of charge bias, the Hamilto- 
nian in the basis of semi-classical states {|Oi), |Oa)> IO3)} 
takes the form Hq — tho r + isthoi with 



iQr 



( cos(nS Q ) cos(nS Q ) 

cos{ttS q ) 1 
y cos(7r,5 Q ) 1 



(46) 



and 



hoi = 



( 

V s * 







(Sf 2 sin(7rf5 Q ) sin(7r<5 Q 
f 2 sin(7r<J Q ) 
)| x sin(7r<5 Q ) 




where 5% = Sf + S] 



The link current operator then takes 
the form (after setting <5* =0; note the sign change for 

the connected tetrahedron, Jj = ~(2e/h)dSHo/d5f) 



lest i 



i(tt5 q ) 




(48) 



In measuring the qubit's state, we first drive the tetra- 
hedron adiabatically towards the measuring point S Q = 
1/2. At small values of 6 = (7r<5 Q ) 2 /2 the doublet splits 
and we find the new eigenvalues Eg and associated eigen- 
vectors \e$) (up to normalization), 



E s /t 



x T 3 

2- # 



a s 



\es) 

•IO2MO3) 
1 ^2 



\ss 



2(l+8/3)\0 1 )-\0 2 )-\0 3 ) 



(49) 



\Q g \ ^ (l-<5/3)|0 1 ) + |02) + |0 3 ) 
\J 3—25/3 



Note that at S — the eigenvectors |a ) and |s ) defining 
the low-energy qubit subspace correspond to the (anti-) 
symmetric combinations |oo) = [|+) — |— )]/v / 2 and |s ) = 
[|+) + |— )]/v / 2- Hence starting with a qubit state |\E') = 
■(/>+!+) ) i n the basis |±), the amplitudes which 

evolve adiabatically towards the measurement point are 
the combinations ip So = + ^_ ) / \[2 and -0 ao = — 
■0-)/\/2 or, in matrix form, 

i>- 
V ^0 

Driving the bias up to 5 Q = 1/2 the spectrum and eigen- 







( ^+ 








V'Oo / 




\ ^0 
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vectors deform into 



-1 \a 1/2 )=i[\0 2 )-\0 3 )]/V2, 

\s 1/2 ) - |Oi>, 

1 |0 va ) = [\0 2 ) + \0 3 )]/V2. 



(50) 



In this new basis {lai/2), l s i/2)i IO1/2)} the current oper- 
ator H48|l takes the form 



2 v / 2est 



/ 1 
1 
yO 



(51) 



hence after the adiabatic evolution, the current mea- 
sures the a x operator in the basis 1(50(1 and has no ma- 
trix elements with the third state |0!/ 2 ) - At the same 
time, the qubit wave function takes the form = 
e-^aok/2} + e w ^ S0 \si/2), where 9 = J dt(E a - 
E s )/h denotes the additional phase picked up in the adi- 
abatic evolution of the qubit amplitudes tp Ss and ^ a6 until 
the measurement is performed. In order to reexpress the 
current 1)51(1 through the original qubit amplitudes ^± 
we define the matrix 



To 



1 

V2 



/ e ~i0/2 _ e -i9/2 q ' 



e i6/2 




je/2 



and find the result 



tfc-i- 



2 v / 2est 



V 




(52) 



describing a high quality measurement of the qubit op- 
erators ±tr z at 9 = 7i7T and ±<r a at 8 = (n + l/2)7r (or 
any component residing in the y-z plane for angles 9 in 
between) . Applying the charge bias 1(44(1 to cyclically per- 
muted islands 1 — ^ 2 — 3 — ^ 1 produces measurements of 
spin-projections in planes rotated by the corresponding 
angles ±27r/3. 

In the following, we restrict the discussion to the mea- 
surement of o z , i.e., we assume that 9 — nir with n even in 
the specific discussion below. Following the result ((5211 . 
the link current J* 1 " 1-1 ' has opposite signs for the two 
qubit-states |+) and |— ) and its measurement allows for 
the determination of the qubit 's final wave function. In 
the actual measurement, the system is driven symmet- 
rically with equal external currents I n i = I m entering 
the system at the points 'nl', 'n2', 'n3' (cf. Fig. 6(a)) 
and leaving symmetrically through the points 'ml', 'm2', 
and 'm3'. The external current I m is chosen close to, 
but below twice the critical current of the large junc- 
tions in the ring, 7 ml < 2/ mc = 4eE m /h. Accounting 
for the additional currents induced in the tetrahedron 
under the charge bias 6 Q = 1/2, the currents /< n2 - ml > 



and 7 <n3 ml) through the large measurement junctions 
are equal to (J m ± 2^/2est/h)/2 for the |±) states of the 



doublet, while the currents I ( 



and I< 



assume 



the values (I m =F 2v / 2esi/fi)/2 for the same eigenstates. 
Note that in the singlet state the bare measurement cur- 
rent / m /2 flows through the junctions as no current is 
induced in the singlet state |0) under the above charge 
bias. As demonstrated in Ref. [3, the measurement cur- 
rent J m and the induced current 2\/2est/h can be chosen 
such that the switching probabilities V± into the tran- 
sient voltage-state of the large measurement junctions are 
strongly different for the |±) states. Then the location 
of voltage pulses on the junctions 'n2-ml' and 'n3-ml' 
or on the complementary junctions l n2-m3' and l n3-m2' 
tells us whether the qubit was in the |+) or |— ) state 
just before the measurement. Furthermore, the absence 
of any voltage pulse is the signature of the singlet state. 

Finally, we discuss the projective measurement of the 
operator a x . Following 1(321) . a flux bias produces finite 
'magnetic fields' h x , y which are bilinear in (5*; one then 
may expect that an appropriate flux bias will produce 
loop currents which are diagonal in the basis \a, s), where 



-))/V2. (53) 



Here, we concentrate on the main contribution (propor- 
tional to v) originating from the modification Vj in the 
energies of the semi-classical minima \Oj). The current 
operators It — —(2e/h)d5Ho/dSf evaluated in the basis 
{]+), |-), |0>} take the form 



2evE 




(54) 



with ei = 4(2<5f + SJ + 5%) and 

Ai = (2£f + £J + <5f)+n/3(5f 

a 3 = _(2a» + a*-a*)-iV3(«? + <5j). 



(55) 



Applying a specific bias 5 2 — ^3 — = ^ produces a 
field along h x (cf. ((32(1 ) and induces the current /< m3 - 3 > = 
h-h = -J< m2 - 2 > in the loop 'ml-m3-3-l-2-m2-ml', 



SevEjS 




(56) 



the current operator j (ml 1} vanishes. Transforming to 
the basis states la, s) of cr T , this takes the form 



r(m3-3) _ 




(57) 



obviously, the anti-symmetric state \a) is already a good 
eigenstate of the current operator, while the states |s) 
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and |0) remain mixed. Diagonalizing, we find the eigen- 
values SevS^j/h and eigenvectors \j), 



J \J) 

1 \a), 

1 [\s)-V2\0)]/V3, 

[V2\s) + |0)]/V3. 



(58) 



The diagonal structure of J (m3 ~ 3) in the subspace spanned 
by | a) and \s) implies a measurement of the spin- 
component a x in the original basis |±). Unfortunately, 
we have not been able to identify a setup producing a 
complete diagonalization of the current operator; still, 
the above constellation allows for a well defined statis- 
tical procedure to identify the qubit state: Let us as- 
sume that just before the measurement the wave func- 
tion of the qubit belongs to the doublet subspace, = 
V> a |a) + ip s \s), with \tp a \ 2 + \4>s\ 2 = 1- The measure- 
ment of the current eigenvalue j then provides us with 
a probabilistic measure of P s = \ip s \ 2 - while the result 
j = appears with a probability (2/3) P g , the current 
value j = 1 is realized with a probability (1/3) P s + P a — 
1 — (2/3) P s , from which the desired quantity P s is eas- 
ily derived within a multi-shot measurement scheme. In 
an experiment tracing the Rabi oscillations between the 
states | a) and \s) the specific time dependence P s (t) = 
cos 2 (f2R, a bii) transforms into the probability difference 
Vj=i(t) ~ Vj=o(t) = 1/3 - (2/3) cos(2f} Rab ii); the reduc- 
tion in the amplitude is a consequence of the admixture 
of the third state |0). 

The actual measurement of the current is carried out 
with the same (symmetric) measurement current config- 
uration as discussed in the previous section. The only 
difference is that in the present measurement the current 
eigenvalues are proportional to j = 1 and j = instead 
of ±1, i.e., the switching of (one of the two) junctions 
'n2-ml', l n3-m2' (or, depending on the sign of flux bias, 
'n2-m3', 'n3-ml') is characteristic for a j = 1 eigenstate, 
whereas the absence of any switching corresponds to the 
j = state. 



V. CONCLUSIONS 

The CH4 molecule can be viewed as a molecular ana- 
logue to the tetrahedral superconducting structure with 
the same tetrahedral symmetry group. The non-Abelian 
character of this symmetry group is responsible for the 
natural appearance of degenerate states. However, con- 
trary to the situation in atomic and molecular physics, 
where such degenerate levels usually correspond to ex- 
cited states, the macroscopic device discussed here can 
be tuned such that the non-Abelian character of its sym- 
metry group manifests itself in the appearance of a de- 
generate ground state. In order to do so, we have to 
bias the device, both electrically and magnetically, to the 
maximally frustrated point with half-flux $0/2 threading 



each sub-loop and half-Cooper-pair charge e induced on 
each island. The ground state is a doublet equivalent to 
a spin-1/2 system in a vanishing magnetic field. 

Given the above analogy to molecular physics, one may 
ask whether a splitting of the doublet ground state due to 
a Jahn- Teller type instability may appear in our system 
as well. In the superconducting tetrahedron this would 
correspond to a paramagnetic instability with a sponta- 
neous breaking of time reversal symmetry: paramagnetic 
currents would lower the system energy due to their in- 
teraction with the self-generated magnetic field and the 
ground state doublet would split — the magnitude of this 
effect is determined by the magnetic inductance of the 
device which we have set to zero in our analysis above. 
However, in order to realize such an instability, the en- 
ergy gain should be linear in the spontaneously generated 
flux; since our device exhibits only a quadratic depen- 
dence on flux, cf. this type of instability is absent. 

On a technical level, the superconducting tetrahedral 
qubit comes with a number of practically useful features: 
i) The weak quadratic sensitivity to electric and magnetic 
noise implies long decoherence times for this device, simi- 
lar to the 'quantronium' discussed by Vion et al£. At the 
same time, the quadratic dependence on bias allows for a 
manipulation via ac-fields rather than the usual dc-bias, 
hence the most dangerous low-frequency part of the noise 
spectrum can be blocked from the qubit via appropriate 
filters, ii) The degenerate ground state allows to avoid 
the appearance of phonon radiation during idle time^: 
Assume a qubit with states |0) and |1) residing at dif- 
ferent energies E§^ E\. A superposition state \^){t) = 
[cxp(-iE o t/h)\0)+f3cxp(-iE 1 t/h)\l)]/ y /l + \P\ 2 will in- 
duce voltage oscillations V — hip/2e cx {E\ — E )/2e 
across the Josephson junction which couple to the un- 
derlying lattice via the piezo-electric effect. Hence, the 
junction acts as an antenna emitting phonons which con- 
tributes to the energy relaxation rate of the qubit. In 
our tetrahedral qubit, E\ = Eq and this loss channel is 
avoided during idle time. Hi) The tetrahedral qubit can 
be fabricated with junctions of relatively large size. This 
is the consequence of a weak cx exp[— const. (Ej/ Ec) 1 ^ 4 ] 
rather than the usual cx exp[— const. (Ej/Ec) 1 ^ 2 ] sup- 
pression of the qubit 's operational energy scale and en- 
tails two important advantages: first, less stringent re- 
quirements on the fabrication process and better junc- 
tion uniformity, and second, an improved robustness of 
the qubit with respect to charge noise originating from 
fluctuating stray charges. The physical origin of this 
benevolent behavior is found in the huge classical ground 
state degeneracy originating from junctions with a simple 
cx sin if current-phase relation combined with a maximal 
magnetic frustration; this degeneracy then is lifted only 
due to quantum fluctuationsi^i^. iv) The tetrahedral 
qubit admits a large variety of manipulation schemes — 
arbitrary manipulations of the effective 'spin 1/2' ground 
state can be implemented through either magnetic or 
electric bias fields, v) The quantum measurement can be 
performed with respect to different basis states and us- 
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ing either charge or flux bias. We have described schemes 
operating in the qubit eigenbasis and projective measure- 
ments onto the current basis. In the former setup, we 
have described detailed procedures for the detection of 
both 'spin' projections a z and a x . In the latter scheme, 
we have identified a high-quality measurement scheme for 
the operators a z and a y through appropriate charge bias; 
the flux bias scheme produces a non-ideal but acceptable 
setup for the measurement of o x . Proper rotation of the 
bias scheme by ±27r/3 provides a measurement of the 
corresponding rotated spin-components. 

The above advantages seem worth the additional com- 
plexity of the device. Still, one may pose the ques- 
tion whether the same benefits can be implemented 
with a simpler device. E.g., the Cs v symmetry group 
of the symmetric three-junction loopS^ also contains a 
two-dimensional representation and appropriate charg- 
ing with qi = 1/3-charge per island produces a doublet 
ground state suitable for quantum computation. How- 
ever, this simpler design does not exhibit the quadratic 
stability under charge noise — charge biasing one island 
reduces the symmetry to Zi with only one-dimensional 
representations and the doublet splits in linear order in 
Sf. On the contrary, charge biasing the tetrahedron re- 



duces the symmetry from to C% v which still contains a 
two-dimensional representation. Indeed, the two complex 
conjugated states |±) react the same way to a charge bias 
Sf and hence the doublet splits only in quadratic order. 
This behavior, the indifferent response of the two ground 
state wave functions to a local perturbation is strongly 
reminiscent of the idea of topological protection, where 
the fault tolerance of the device is implemented on the 
hardware— rather than the software level^. The above 
symmetry arguments then show, that in order to ben- 
efit from a protected degenerate ground state doublet, 
the qubit design requires a certain minimal complexity; 
it seems to us that the tetrahedron exhibits the minimal 
symmetry requirements necessary for this type of protec- 
tion and thus the minimal complexity necessary for its 
implementation. 
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